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Abstract 

We study the radiative 7r~^7r~^ and decays and we calculate 

their branching ratios using a phenomenological approach by adding to the 
amplitude calculated within the framework of chiral perturbation theory and 
vector meson dominance the amplitude of fi-meson intermediate state. Our 
results for the branching ratios are in good agreement with the experimental 
values. 
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The radiative decays of neutral vector mesons into a single photon and a pair of neutral 
pseudoscalar mesons have been a subject of continuous interest. The studies of such decays 
may serve as tests for the theoretical ideas about the nature of the intermediate states and 
the interesting mechanisms of these decays, and they may thus provide information about 
the complicated dynamics of meson physics in the low energy region. 

The very recent measurement of the branching ratio for the decay p'^ tt'^tt'^^ by the 
SND Collaboration obtained the value BR{p'^ 7r°7r%) = (4.llJ:^ ± 0.3) x lO^^ [1], thus 
improving their previous preliminary report of BR{p^ 7r°7r°7) = (4.8li;8 i 0-3) x 10~^ 
[2]. On the other hand, the branching ratio for the decay p° — > 7r"'"7r~7 was reported earlier 
by Novosibirsk group as BR{p^ 7r"'"7r~7) = (9.9 ± 1.6) x 10~^ [3,4], and it was observed 
that the pion bremsstrahlung is the main mechanism for this decay with the structural 
radiation proceeding through the intermediate scalar resonance making less than one-order 
of magnitude smaller contribution to the branching ratio [3]. 

The theoretical studies of radiative p-meson decays was initiated by Singer [5] who cal- 
culated the amplitude for the decay p° — > 7r+7r~7 by considering only the bremsstrahlung 
mechanism, and he assumed that the decay p° — > 7r°7r°7 proceeds through an (a;7r) interme- 
diate state as p^ — > {uj)t^^ — >■ (7r°7)7r'^. The vector meson dominance (VMD) calculation of 
Bramon et al. [6] with this intermediate state using standard Lagrangians obeying SU(3)- 
symmetry gave the value B(p° 7r°7r°7)=l.l x 10~^ for the branching ratio. However, 
they also noted that final state interactions could lead to a larger value for the branch- 
ing ratio B(p° ir^n^'-f) through the mechanism p° — >■ (7r"'"7r~)7 — * (71^11^)'^. Bramon et 
al. [7] later considered the radiative vector meson decays within the framework of chiral 
effective Lagrangians enlarged to include on-shell vector mesons using chiral perturbation 
theory, and they calculated the branching ratios for various decays at the one-loop level, 
including both tttt and KK intermediate loops. In this approach, the decay p° 7r°7r°7 
proceeds mainly through the charged pion (7r''"7r'~)loops, contribution of kaon- loops being 
three orders of magnitude smaller, resulting in the decay rate r(p'^ 7r°7r'^7)^=1.42 KeV 
which is of the same order of magnitude as the VMD contribution. The interference between 
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the pion-loop contribution and the VMD-amphtude turns out to be constructive leading to 
BR(p° 7r°7r°7)yMD+x=2-6 x 10~^. However, this value is still substantially smaller than 
the latest experimental result quoted above. 

Since the experimental result for the branching ratio B(p° — > 7r°7r°7) is almost nearly 
twice the theoretical value calculated using VMD and chiral-loop amplitudes, the mechanism 
of the decay p° tt^tt^j should reexamined and additional contributions should be investi- 
gated. One additional contribution to the decay may be provided by the amplitude involving 
scalar-isoscalar cr-meson as an intermediate state. Although the existence of cr-meson has 
long been controversial, an increasing number of theoretical and experimental analyzes find a 
(T-pole position near (500 — i250) MeV [8]. Furthermore, the — > (77r° — > Stt decay channel 
observed by the Fermilab (E791) Collaboration is interpreted to provide direct experimental 
evidence for cr-meson where it is seen as a clear dominant peak with M^. = (478l|3 =t 17) 
MeV, and Fo- = (324l4o ± 21) MeV [9]. Since cr-meson is assumed to couple strongly to 
low mass pion pairs, the p° — > tttt^ decays thus provide an opportunity to investigate the 
theoretical ideas about the role cr-meson plays in the dynamics of low energy meson physics. 

One way to include the effect of cr-meson in the decay mechanisms of p° — > tt'^tt''^ and 
p° — > 7r°7r°7 decays is to consider its contribution as resulting from a cr-pole intermediate 
state, that is to assume that the contributions of cr-meson to the decay mechanisms of 
these decays result from the corresponding amplitudes of {aj) (7r"'"7r~)7 and p° — > 
(0-7) — > (7r°7r°)7 reactions. In a previous work [10], two of the present authors calculated the 
branching ratio B{p^ — >■ 7r"'"7r~7)in a phenomenological framework using pion bremsstrahlung 
amplitude and cr-meson pole amplitude. The experimental value of this branching ratio was 
then used to calculate the coupling constant gp^-y as a function of cr-meson parameters 
and Fo-. These authors in a following work [11] calculated the branching ratio B(p° — > 7r°7r°7) 
using the values of the coupling constant gp^j thus obtained in a phenomenological approach 
where the contributions of cr-meson, w-meson intermediate states and of the pion-loops are 
considered. The branching ratio B{p^ — >■ 7r''7r''7) obtained this way for M^j = 478 MeV 
and Fo- = 324 MeV was more than an order of magnitude larger than the experimental 
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value. This unrealistic value was the result of the constant p ^ amplitude used and 
consequently the large coupling constant Qp^^ extracted from the experimental value of the 
branching ratio of the — > 7r+7r~7 decay. Therefore, it may be concluded that it is not 
realistic to include the cr-meson in the mechanisms of the radiative p'^-meson decays as an 
intermediate pole state. 

On the other hand, p° — > 7r°7r°7 decay was also considered by Marco et al. [12] in 
the framework of unitarized chiral perturbation theory. They noted that the energies of 
two pion system are quite large so that the decay cannot be treated with standard chiral 
perturbation theory. They used the techniques of chiral unitary theory developed earlier, 
a review is given by Oiler et al. [13], to include the final state interactions of two pions 
by summing the pion-loops through the Bethe-Salpeter equation. The branching ratio for 
p° — >• 7r°7r°7 they obtained was B(p° — >• 7r°7r°7) = 1.4 x 10~^ and they furthermore noted that 
this result could be interpreted as the result of the mechanism — > ((7)7 — > (7r°7r°)7 since 
7r°7r° interaction is dominated by cr-pole in the relevant energy regime of this decay. 

Thus, it seems that a natural way to include the effects of cr-meson in the mechanisms 
of radiative p°-meson decays is to assume that the cr-meson couples to p°-meson through 
the pion-loop. In this work, we reconsider the approach used in the references [10] and [11], 
and we study the contribution of the cr-meson intermediate state amplitude to — > 7r"'"7r~7 
and p° t^^t:^^ decays. We follow a phenomenological approach and assume that the 
cr-meson couples to the p°-meson through a pion-loop, in other words we assume that the 
amplitude p° — >• 0^7 results from the sequential p° {ti^ti~)'^ — >• cr7 mechanism as suggested 
by the unitarized chiral perturbation theory in which cr-meson is generated dynamically by 
unitarizing the one-loop pion amplitudes. We use the coupling constants that are determined 
from the experimental values of the relevant quantities that we calculate employing the 
effective Lagrangians of our approach. Although the decay p° — > 7r"'"7r~7 is dominated by 
the pion-bremsstrahlung amplitude, for the consistency of our formalism we study both of 
p° — > 7r"'"7r~7 and p" tt^tt^'j decays and calculate their branching ratios. 

Our calculation is based on the Feynman diagrams shown in Fig. 1 for p° 7r"'"7r~7 
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decay and in Fig. 2 for p° — > 7r°7r°7 decay. The last diagrams in Fig. 1 a, b, c and in 
Fig. 2 b, c are the direct terms required to estabhsh the gauge invariance. We describe the 
a; pTT- vertex by the effective Lagrangian 

which also conventionally defines the coupling constant guipn- Achasov et al. [14] assumed 
that a; — > Stt decay proceeds with the intermediate pir state as a; — > (p)7r — > tttttt and 
using experimental value of the cu ^ Stt width they determined this coupling constant as 
gpuiiz = (14.3 ± 0.2) GeV~^. The ptttt- vertex is described by the effective Lagrangian 

C%i = g,^^p,-{d^7rX7r) . (2) 

The experimental decay width of the decay p — > tttt [4] then yields the value gip^Tr = (6.03 ± 
0.02) for the coupling constant gf^n- For the crTTTr-vertex we use the effective Lagrangian 

Clii ^ ^g^^^M^n ■ na . (3) 

Using the experimental values for and [9] as = (483±31) MeV and F^ = (338±48) 
MeV, where statistical and systematic errors are added in quadrature [15], we obtain the 
strong coupling constant g^T^T: as gann — (5.34±0.55). We note that our effective Lagrangians 
£^{4 and are the ones that result from an extension of the cr-model to include the 
isovector p through a Yang-Mills local gauge theory based on isospin with the vector meson 
mass generated through the Higgs mechanism [16]. The a;7r7- vertex is described by the 
effective Lagrangian 

Ct/J^^g^^^e'''"'fd^u;,d^A0n . (4) 

We then obtain the coupling constant g^^T^^ from the experimental partial width [4] of the 
radiative decay lo 7r°7 as gu^-n-y — (0.706 ± 0.021) GeV~^. 

Meson-meson interactions were studied by Oiler and Oset [17] using the standard Chiral 
Lagrangian in lowest order of chiral perturbation theory that contains the most general low 
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energy interactions of the pseudoscalar meson octet in this order. We use their results for 
the four pseudoscalar amplitudes tt+tt" 7r°7r° and tt+tt" tt'^tt^ that we need in the 
loop diagrams in Fig. 1 b and in Fig. 2 b. We note that as shown by Oiler [18] due 
to gauge invariance the off-shell parts of the amplitudes, that should be kept inside the 
loop integration, do not contribute, and as a result the amplitudes AlT^(7r'^7r~ — > 7r°7r°) and 
^^^(Tr+Tr" TT+TT") factorize in the expressions for the loop diagrams. 

In our calculation of the invariant amplitude, in the cr-meson propagator we make the 



replacement — 
which is given as 



g,2 _ M^ + iMT and we use the energy dependent width for cr-meson 



^a{q ) = 



A,' -AMI) 



(5) 



\Ml- AMI 

For the loop integrals appearing in Figs. 1 and 2 we use the results of Lucio and Pestiau 
[19] who evaluated similar integrals using dimensional regularization. In our case the con- 
tribution of the pion-loop amplitude corresponding to p° (7r^7r^)7 — > 7r°7r°7 reaction can 
be written as 



= 27rW^ ■k){e-u)-{p- e){k ■ u)] , 



(6) 



where a = M^/M^, b = {p - kf/M^, M^{7r+7r- tt^tt^) 



is - MD/fl s = 



M^o^o, /tt = 92.4 MeV, p{u) and k{e) being the momentum (polarization vector) of p- 
meson and photon, respectively. The amplitudes corresponding to p° (7r"'"7r~)7 — > 7r'''7r~7, 
p° — >• (7r+7r~)7(T — >• 77r+7r~ and p° — >• (7r"'"7r~)7cr — >• 77r°7r° reactions can similarly be written. 



The function I(a,b) is given as 
1 2 



7(a,6) = 



2(a -h) (a - hf 



+ 



[a - by 



(7) 



where 



- [arcsm^^v: 



X > 
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{Ax - 1)2 arcsin(2^) , x > j 



i(l-4a:)i[ln(^)-Z7r 

l±(l-4x)5] . (8) 



We then calculate the invariant amphtude A1(E^, Ei) from the corresponding Feynman di- 
agrams shown in Fig. 1 and 2 for the decays p° — > 7r+7r~7 and p° — > 7r°7r°7, respectively. The 
differential decay probability for an unpolarized p°-meson at rest is then given as 

= M (9) 

dE^dEi (27r)3 8Mp ' ' ' ^ ^ 

where and Ei are the photon and pion energies respectively. We perform an average over 
the spin states of p°-meson and a sum over the polarization states of the photon. The decay 
width is then obtained by integration 

.max. (] I 

/ dE^ dEi " (10) 

,Tnin. (JiJ-J^UjJ-J^ 

where the factor (|) is included for the calculation of the decay rate for p° — > 7r°7r°7 because 
of the TT^TT*^ in the final state. The minimum photon energy is E^ „^j„. = and the maxi- 
mum photon energy is given as Ej^j^ax. — (-^p — 4M^)/2Mp=338 MeV. The maximum and 
minimum values for pion energy Ei are given by 

^ -[-2E^Mp + 3E^M^ -M^ 



±E^^{-2E^M, + Ml){-2E^Mp + - AMI) 



The photon spectra for the branching ratio of the decay p° tt^tt'j are plotted in Fig. 
3 as a function of photon energy E^. The contributions of bremsstrahlung and structural 
radiation amplitudes calculated with pion-loop and with cx-meson intermediate state as well 
as the contribution of the interference term are shown as a function of the photon energy, 
where the minimum photon energy is taken as E^^min — 50 MeV since the experimental 
value of the branching ratio is determined for this range of photon energies [3] . On the same 
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figure we also show the experimental data points taken from Ref. [3] which are normalized 
to our results. As it can be seen in Fig. 3, the shape of the photon energy distribution is 
reproduced well. As expected, the main contribution to the branching ratio comes from the 
pion-bremsstrahlung amplitude, the contributions of pion-loop and cr-meson intermediate 
states becoming noticeable only in the region of high photon energies. On the other hand, 
if a cr-meson pole model or equivalently a constant p — > (77 amplitude is used as in Ref. [10] 
the contribution of the u-term becomes increasingly important in the region of high photon 
energies dominating the contribution of the bremsstrahlung amplitude, and although its 
contribution is somewhat reduced by the interference term the cr-meson amplitude makes 
the main contribution to the branching ratio in the region of high photon energies confiict- 
ing with the experimental spectrum. For the contribution of different amplitudes to the 
branching ratio, we obtain B(p° — >■ 7r"'"7r~7)-y = (1.14 ±0.01) x 10~^ from the bremsstrahlung 
amplitude, B(p° 77+77^7)^ = (0.45 ± 0.08) x 10~^ from the pion-loop amplitude and 
B(p° — > 7r+7r~7)(^ = (0.83 ± 0.16) x 10~^ from the u-meson intermediate state. If we con- 
sider also the interference between the 7r-loop and cr-meson amplitudes, we then obtain the 
contribution coming from the structural radiation as B(p° — > 7r"'"7r~7) = (0.83 ±0.14) x 10~^ 
which is reasonable agreement with the experimental limit B(p'^ — > 7r'*'7r~7) < 5 x 10~^ 
deduced by Dolinsky et al. [3] for the structural radiation. Dolinsky et al. [3] also extracted 
the experimental limit B(p° e(700)7 n~^7r~'-f) < 4 x 10~^ where the transition proceeds 
through the intermediate scalar resonance which is considerably lower. Our result for the 
contribution of the cr-meson intermediate state B(p° n^n^'y)^ = (0.83 ± 0.16) x 10~^ is 
also in good agreement with this experimental limit. For the total contribution, we obtain 
the branching ratio of the decay p° 7r+7r~7 as B(p° 7r+7r~7) = (1.22 ± 0.02) x 10~^ for 
Ey > 50 MeV. Our result is in reasonably good agreement with the experimental number 
B(p° ^ 7r+7r-7) = (0.99 ± 0.16) x 10"^ [3]. 

The photon spectra resulting from our calculation for the branching ratio of the p° 
7r°7r°7 decay is shown in Fig. 4. The contributions of VMD-amplitude, pion-loop amplitude 
and cr-meson intermediate state amplitude as well as the total interference term are indicated. 
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The contributions of VMD amplitude, pion-loop amplitude and cr-meson intermediate state 
amplitude to the branching ratio of the decay are B(p° '7^^'7r°j)vMD = (1.03±0.02) x 10~^, 
B(pO ^ 7rV7)^ = (1.07 ± 0.02) x lO'^ and B(pO ^ 7r°7r°7)<, = (4.96 ± 0.18) x 10"^ 
respectively. We see that cr-meson intermediate state makes an important contribution to 
the branching ratio comparable to the contributions of VMD- and pion-loop amplitudes. 
The values we obtain for B(p° —>■ 7r°7r°7)yMD and for B(p° — > 7r°7r°7)7r are in agreement 
with previous calculations [6,7]. For the total branching ratio, including interference terms, 
we obtain the result B(p° — >^ 7r°7r°7) = (4.95 ±0.82) x 10~^ which is in a better agreement 
with the experimental result B{p^ — > tt^tt °7) = (4.11^:^ ±0.3) X 10-5 [1] than the theoretical 
value obtained by using only VMD- and chiral pion-loop amplitudes [7]. 

In our work, we follow a phenomenological approach and in a model for the decay mech- 
anism of p°-meson radiative decays including the contribution coming from the cr-meson 
intermediate state as well as VMD- and chiral pion-loop contributions we calculate the 
branching ratios for the p° — > n'^n'^ and p° 7r°7r°7 decays. In our calculations of the 
branching ratios, the coupling constants that we use in our model are determined from the 
relevant experimental quantities. Our results for the branching ratios are in good agree- 
ment with the experimental values, and we believe that our study demonstrates that the 
contribution coming from the cr-meson intermediate state amplitude should be included in 
the analysis of radiative p°-meson decays, and moreover cr-meson should be considered to 
couple to the p° meson through a pion-loop. 

In a recent paper, Palomar et al. [20] also evaluated the branching ratios of the radiative 
p° and cu decays into 7r°7r° and tt^t]. They used the sequential vector decay mechanisms 
in addition to chiral loops and p — uj mixing. For the p° — > 7r°7r°7 decay p — uj mixing 
was negligible, but the branching ratio obtained with the sum of the sequential and loop 
mechanisms was about three times larger than with either mechanism alone, leading to a 
result -B(p° n^n^j) — 4.2 x 10"^ comparable with the present experimental value. 

Therefore, in order to understand the mechanism of the p° — >■ 7r°7r°7 decay and to obtain 
insight into the nature and the properties of the cr-meson, and the role it plays in the 
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dynamics of low energy meson physics, further experimental tests such as the measurements 
of invariant mass distributions will be very valuable. 
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FIGURES 




FIG. 1. Feynman diagrams for the decay p n'^n 7. 
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FIG. 2. Feynman diagrams for the decay 7r°7r°7. 
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FIG. 3. The photon spectra for the branching ratio of 7r^7r~7 decay. The contributions 

of different terms are indicated. The experimental data taken from Ref. [3] are normaHzed to our 
results. 
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FIG. 4. The photon spectra for the branchmg ratio of decay. The contributions of 

different terms are indicated. 
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